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A SPECTRAL METHOD FOR CALIBRATING A MULTI-AXIS 
ACCELEROMETER DEVICE 

ABSTRACT OF THE DISCLOSURE 

A method for simultaneously determining respective scale factors and 
alignment angles of a multi-axis accelerometer device for measuring 
acceleration. To measure the scale factors, the multi-axis accelerometer 
device to be calibrated is mounted on a turntable. The turntable has a tilt 
angle with respect to a vertical axis defined by the local gravity vector. The 
turntable is used to spin the multi-axis accelerometer device around an axis of 
rotation at an angular velocity such that the multiple sensitive axes of the 
accelerometer device experience a time varying component of the local gravity 
vector. The respective outputs of the multiple sensitive axes of the 
accelerometer device are recorded as the device experiences the time varying 
component of the local gravity vector. The multi-axis accelerometer device is 
mounted on the turntable in two more orthogonal orientations, and the rotation 
and data recording procedure is repeated. The angular velocity can be constant 
during the logging. Additionally, the predicted output of an ideal accelerometer 
on the turntable is generated, wherein the predicted output is based on the tilt 
angle of the turntable and the angular velocity of the turntable and on the 
value of gravitational acceleration at the location of calibration. The data 
recorded with the accelerometer device mounted in the three orientations are 
combined with the predicted output, to obtain the scale factors and alignment 
angles of the multi-axis accelerometer device. 
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This application is related to copending application "A SPECTRAL METHOD 
FOR CALIBRATING ACCELEROMETERS" by Nadkarni et al., filed on 
( serial No. , which is incorporated herein by reference. 



A RPFP.T RAT, METHOD FOR CALIBRATING A MULTI-AXIS 
AfiflF.T UROMETER DEVICE 



10 TECHNICAL FIELD 

The present invention relates generally to methods and systems for the 
5 precise calibration of instruments. More specifically, the present invention 
S pertains to an accurate and efficient process for calibrating accelerometers. 
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gjl5 BACKGROUND ART 

An accelerometer is a transducer used for measuring acceleration. 
Acceleration is usually measured at a measurement point in the 
W accelerometer, along a sensitive axis of the accelerometer. Generally, the 

magnitude of an applied acceleration is communicatively coupled to external 
instruments or circuits as an electrical impulse having amplitude proportional 
to the magnitude of the applied acceleration. The electrical impulse comprises 
the measured acceleration and is processed by the external circuits as required 
for a variety of applications. One such application is, for example, an Inertial 
Measurement Unit (IMU), where acceleration measurements are used to 
25 generate velocity and positioning information. 

The electrical impulse output of an accelerometer is proportional to the 
acceleration, applied at the measurement point along the sensitive axis of the 
accelerometer. The process of calibrating an accelerometer consists of 
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computing a constant of proportionality, referred to as a scale factor of the 
accelerometer. The scale factor of an accelerometer precisely relates the 
amplitude of the electrical impulses comprising the measured acceleration to 
the magnitude of a corresponding acceleration applied at the measurement 
point, along the sensitive axis of the accelerometer. 

A multi-axis accelerometer device can measure acceleration along 
multiple sensitive axes. This can be a combination of one or more 
accelerometers, with one or more axes of sensitivity each, and a common 
frame of reference with respect to which each of these accelerometers and 
their respective measurement points and sensitive axes remains fixed at all 
times. The frame of reference of the multi-axis accelerometer device is the 
coordinate system in which the acceleration sensed by the array is measured. 
The frame of reference of the array is typically an orthogonal frame of 
reference. 

The process of calibrating a multi-axis accelerometer device consists of 
computing the scale factors for each of the multiple sensitive axes in the 
device, and furthermore computing the alignment angles of the sensitive axes 
in the device with respect to a frame of reference of the device. One measure of 
the alignment angles of a sensitive axis of the device is the direction cosine 
vector or alignment vector of this sensitive axis of the accelerometer device 
with respect to the orthogonal frame of reference of the array. The alignment 
vector of a sensitive axis of the multi-axis accelerometer device is the unit 
vector in the direction of the sensitive axis of the device. For optimal precision 
of measurement using the multi-axis accelerometer device, it is desirable to 
calibrate the multi-axis accelerometer device by precisely determining the 
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scale factors and alignment angles corresponding to each individual sensitive 
axis of the device. 



Prior art systems for calibrating accelerometers (e.g., measuring and 
defining the scale factor) relied on comparisons of the accelerometers to certain 
10 standard devices. Such prior art systems necessarily assume that the 
standard devices themselves are properly calibrated, often leading to the 
introduction of additional error into the calibration process. For example, one 
prior art system (see prior art US Patent 5,970,779) requires the use of 

0 precisely controlled swing arm motor systems to which the accelerometer 

1 15 being tested is mounted, along with an appropriate counter weight. The swing 
arm motor would be precisely controlled by a processor to impart a simple 
harmonic motion acceleration to the sensitive axis of the accelerometer, and 
vary this acceleration by varying the angular acceleration of the swing arm. 
The resulting output of the accelerometer would be examined with respect to 

20 the controlled varying of the swing arm motor, and the scale factor would be 
determined therefrom. 
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One problem with the above prior art approach is that it requires a 
precisely controllable motor for varying the angular velocity of the 
25 accelerometer. The motor needs to precisely apply a simple harmonic 

acceleration to the accelerometer by varying the angular velocity about an 
axis of rotation. As described above, this system requires the proper 
calibration of the standard devices themselves (e.g., the motor), which often 
leads to additional error in the calibration of the accelerometer. 

30 
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A second, more important drawback of the above prior art approach is 
that it requires measuring the radius of rotation of the accelerometer. This 
distance can be very difficult to measure accurately, since the measurement 
point of the accelerometer is internal to the accelerometer. Any error in this 
measurement will manifest itself in through a flawed calibration. 

Thus, what is required is a solution that accurately measures and 
determines the scale factor and alignment angles of each of the multiple 
sensitive axes in the device simultaneously. What is required is a solution that 
calibrates the multi-axis accelerometer device without introducing 
unnecessary sources of error. The required solution should be precise and avoid 
reliance on standard devices, which can introduce error into the calibration 
process. The required solution should not rely on any time varying control of a 
standard device to impart variable acceleration. The required solution should 
not rely on measurements of distance to points internal to the accelerometer. 
The present invention provides a novel solution to the above requirements. 
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5 DISCLOSURE OF THE INVENTION 

The present invention provides a solution that accurately measures and 
determines the scale factor and alignment angles of multiple sensitive axes of a 
multi-axis accelerometer device simultaneously. The present invention 
provides a solution that calibrates a multi-axis accelerometer device without 
10 introducing unnecessary sources of error. The solution of the present invention 
is precise and avoids reliance on standard devices, which can introduce error 
into the calibration process. The solution of the present invention does not rely 
Q on any time varying control of a standard device to impart variable 



N 15 



acceleration. 

In one embodiment, the present invention is implemented as a spectral 
method for simultaneously determining respective scale factors and alignment 
vectors of a multi-axis accelerometer device for measuring acceleration. The 
scale factors and alignment angles are determined simultaneously in one 



aft « 

M 

^ 20 process, allowing the calibration of the multiple axes of the multi-axis 



accelerometer device in one process. To measure the scale factors and 
alignment angles, the multi-axis accelerometer device to be calibrated is 
mounted on a turntable. The turntable has a tilt angle with respect to a 
vertical axis defined by the local gravity vector. The turntable is used to spin 

25 the multi-axis accelerometer device around an axis of rotation at an angular 
velocity such that each sensitive axis of the device experiences a time varying 
component of the local gravity vector (e.g., due to the tilt angle). The 
respective outputs of the multiple sensitive axes of the multi-axis 
accelerometer device are logged as each sensitive axis in the array experiences 

30 the time varying component of the local gravity vector. This process is 
repeated with the multi-axis accelerometer device placed in each of three 
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orthogonal orientations along the axes of the frame of reference of the array 
(e.g. the orthogonal X, Y and Z axes of the frame of reference). 

The scale factors and alignment vectors of the sensitive axes of the 
multi-axis accelerometer device are determined by combining the recorded 
outputs of the multiple sensitive axes of the device mathematically with the a 
predicted output of an ideal accelerometer (e.g., a sine wave). Herein, the 
predicted output is based on the tilt angle, the angular velocity of the ideal 
accelerometer and on gravitational acceleration. This combination is 
performed after the multi-axis accelerometer device has been placed on the 
turntable in each of three orthogonal orientations along the axes of the frame 
of reference of the device. In so doing, the present invention accurately 
measures and determines the scale factor and alignment angles of the 
sensitive axes of the device, without relying on any time varying control of a 
standard device (e.g., stepper motor, etc.) to impart variable acceleration to 
the multi-axis accelerometer device. It also does not rely on a measurement of 
the radius about which the multi-axis accelerometer device rotates, and 
thereby, does not rely on a precise knowledge of the location of measurement 
point of the device internal to the device. 
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5 BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying drawings, which are incorporated in and form a part 
of this specification, illustrate embodiments of the invention and, together with 
the description, serve to explain the principles of the invention: 

10 Figure 1A shows a multi-axis accelerometer device calibration system in 

accordance with one embodiment of the present invention. 

Figure IB shows a diagram of 3 sensitive axes (e.g. A, B and C) of the 
multi-axis accelerometer device from Figure 1A, along with the frame of 
reference (e.g. axes X, Y and Z) of the multi-axis accelerometer device. 

Figure 1C shows a diagram of the device in Orientation 1, in which the Z- 
axis of the frame of reference of the device is pointing along the axis of rotation 
of the turntable. 

Figure ID shows a diagram of the device in Orientation 1, in which the Y- 
axis of the frame of reference of the device is pointing along the axis of rotation 
of the turntable. 

25 Figure IE shows a diagram of the device in Orientation 1, in which the X- 

axis of the frame of reference of the device is pointing along the axis of rotation 
of the turntable. 



Figure 2 shows the components of the multifunction processor in 
30 accordance with one embodiment of the present invention. 
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Figure 3 shows a flowchart of the steps of an accelerometer calibration 
process in accordance with one embodiment of the present invention. 

Figure 4 shows a diagram of the basic components of a computer 
system in accordance with one embodiment of the present invention. 



TRMB 964 



8 



July 27, 2001 



5 DETAILED DESCRIPTION OF THE INVENTION 

Reference will now be made in detail to the embodiments of the 
invention, a spectral method for calibrating a multi-axis accelerometer device, 
examples of which are illustrated in the accompanying drawings. While the 
invention will be described in conjunction with the preferred embodiments, it 
10 will be understood that they are not intended to limit the invention to these 

embodiments. On the contrary, the invention is intended to cover alternatives, 
modifications and equivalents, which may be included within the spirit and 
scope of the invention as defined by the appended claims. Furthermore, in the 

0 

£ following detailed description of the present invention, numerous specific details 
£ 15 are set forth in order to provide a thorough understanding of the present 
'i invention. However, it will be obvious to one of ordinary skill in the art that the 
tj! present invention may be practiced without these specific details. In other 



instances, well-known methods, procedures, components, and circuits have not 
been described in detail as not to unnecessarily obscure aspects of the present 



determines the scale factor and alignment angles of each of the multiple 
sensitive axes of the accelerometer device. The present invention provides a 

25 solution that can calibrate a device having multiple sensitive axes in a single 
calibration process. The present invention provides a solution that calibrates 
accelerometer devices without introducing unnecessary sources of error. The 
solution of the present invention is precise and avoids reliance on standard 
devices, which can introduce error into the calibration process. The solution of 

30 the present invention does not rely on any time varying control of a standard 
device to impart variable acceleration. The solution of the present invention 




invention. 



The present invention provides a solution that accurately measures and 
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5 does not rely on any measurement of distance to a measurement point of the 
device. The present invention and its benefits are further described below. 

Notation and Nomenclature 

Some portions of the detailed descriptions, which follow, are presented in 
10 terms of procedures, steps, logic blocks, processing, and other symbolic 

representations of operations on data bits within a computer memory. These 
descriptions and representations are the means used by those skilled in the 
data processing arts to convey most effectively the substance of their work to 
0 others skilled in the art. A procedure, computer-executed step, logic block, 
W15 process, etc., are here, and generally, conceived to be self-consistent sequences 
H of steps or instructions leading to a desired result. The steps are those 
W requiring physical manipulations of physical quantities. Usually, though not 
s necessarily, these quantities take the form of electrical or magnetic signals 

vjj capable of being stored, transferred, combined, compared, and otherwise 

Piii 

ij|20 manipulated in a computer system. It has proven convenient at times, 
v principally for reasons of common usage, to refer to these signals as bits, 
values, elements, symbols, characters, terms, numbers, or the like. 

It should be borne in mind, however, that all of these and similar terms 
25 are to be associated with the appropriate physical quantities and are merely 

convenient labels applied to these quantities. Unless specifically stated 

otherwise, as apparent from the following discussions, it is appreciated that 

throughout the present invention, discussions utilizing terms such as 

"processing," "computing," "configuring," "comparing," "detennining," 
30 "sampling," "transforming," or the like, refer to the action and processes of a 

computer system (e.g., computer system 404 of Figure 4), or similar electronic 
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computing device, that manipulates and transforms data represented as 
physical (electronic) quantities within the computer system's registers and 
memories into other data similarly represented as physical quantities within 
the computer system registers or memories or other such information storage, 
transmission, or display devices. 

Embodiment of the Invention 

Referring now to Figure 1A, a multi-axis accelerometer device 
calibration system 100 in accordance with one embodiment of the present 
invention is shown. As depicted in Figure 1A, the calibration system 100 
utilizes a turntable 130 and a multifunction processor system 120 to calibrate 
the output of a multi-axis accelerometer device 110. The multi-axis 
accelerometer device 110 is mounted on turntable 130 in one of three 
orientations, shown in Figure 1C, ID and IE, and described below. The axis of 
rotation of turntable 130 is tilted with respect to a vertical Z axis by an angle 
9. The Z axis is precisely vertical, parallel to the local gravity vector, g. During 
the calibration process, the turntable rotates multi-axis accelerometer device 
110 about the axis of rotation at an angular velocity co. It should be noted that 
the Z axis depicted in Figure 1A refers to the vertical axis, parallel to the local 
gravity vector, as opposed to any Z axis of the multi-axis accelerometer device 
110. 

The rotation of the turntable about the axis of rotation at the constant 
angular velocity co gives rise to three forces acting on the multi-axis 
accelerometer device 110, shown as Fl, F2 and F3. F2 is a centripetal force of 
constant magnitude acting on the multi-axis accelerometer device. The 
gravitational acceleration, not shown in the figure, felt by the multi-axis 
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5 accelerometer device in the plane of the turntable is g*sin(9). Decomposing this 
acceleration along the radius of rotation and tangential to the radius of rotation 
gives us forces Fl and F3. Fl is a time varying force with amplitude 
g*sin(9)*cos((|)(t)) acting radially on the multi-axis accelerometer device, and F3 
is a time varying force with amplitude g*sin(9)*sin(<|)(t)), acting tangentially on 
10 the multi-axis accelerometer device. Here 9 is the angle of tilt of the axis of 
rotation from the vertical (e.g., the Z axis), as shown in Figure 1A. <|) is the 
angle subtended at the axis of rotation by the frame of reference of the multi- 
axis accelerometer device 110 and the component of gravity in the plane of 
? rotation of the accelerometer 110. <|> is a function of time t. g is the 

f? 15 acceleration due to gravity. The angular velocity co=d<t>/dt. 

M ! 

H Specifically, <|)(t) = co*t + <[>(0), where <|>(0) is the value of the angle § at time t=0, 

W which is when data logging begins. 

W 

\| Referring now to Figure IB, the multi-axis accelerometer device in this 

?i ij 

\]i 20 embodiment consists of three accelerometers A, B and C, which are oriented in 

o 

a fixed orientation with respect to each other. Typically, these accelerometers 
A, B and C are nominally mutually orthogonal. Any mounting error present is 
compensated for when the alignment angles of the accelerometer are 
computed. Additionally, the multi-axis accelerometer device has an orthogonal 

25 frame of reference with respect to which it measures acceleration. Each of the 
accelerometers A, B and C are also oriented in a fixed orientation with respect 
to this orthogonal frame of reference. Typically, the axes of the orthogonal 
frame of reference of the multi-axis accelerometer device (e.g. X, Y and Z axes 
of the device) nominally coincide with the sensitive axes of the three 

30 accelerometers. Thus, accelerometer A is aligned to point in the X direction of 
the multi-axis accelerometer device, but is not perfectly aligned. Similarly, 
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5 accelerometer B and C are aligned, but not perfectly, to respectively point 
along the Y and Z axes of the orthogonal frame of reference of the multi-axis 
accelerometer device 110. It is necessary to determine the alignment angles of 
the accelerometers, because the accelerometers A, B and C are not perfectly 
aligned with the orthogonal frame of reference of the multi-axis accelerometer 
10 device. 



Referring again to Figure IB, the multi-axis accelerometer device 110 is 
shown with its multiple sensitive axes (e.g., along the arrows A, B and C) which 
are the sensitive axes of the individual accelerometers (e.g., accelerometers A, 
15 B and C) in the device. As shown in Figure IB, the multi-axis accelerometer 
J device 110 includes three sensitive axes of measurement. The multiple axes 

3 . ; : 

fi are nominally mutually orthogonal, and are used to provide acceleration 

® measurements, and position derived therefrom, in 3D space. Each axis has its 

Q 

ht j corresponding output, shown as V B and V c . In the present embodiment, the 

m 

* ff | 20 outputs convey the measured acceleration as voltage signals. A connector CI 
M> couples these output voltages V A V B and V c from the multi-axis accelerometer 
device 110 to the multifunction processor 120 as depicted in Figure 1A. In this 
embodiment, a connector C2 connects the multifunction processor system 120 
to turntable 130 to receive information regarding the angular velocity co. 

25 

Referring now to Figures 1C, ID and IE, these figures show the multi- 
axis accelerometer device positioned in each of three orientations. In 
Orientation 1, shown in Figure 1C, the Z-axis of the multi-axis accelerometer 
device is parallel to the Axis of Rotation of the turntable, and the X-axis of the 
30 multi-axis accelerometer device is aligned radially, pointing along the direction 
of force Fl. In Orientation 2, shown in Figure ID, the Y-axis of the multi-axis 
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accelerometer device is parallel to the Axis of Rotation of the turntable, and 
the X-axis of the multi-axis accelerometer device is aligned radially, pointing 
along the direction of force Fl. In Orientation 3, shown in Figure IE, the X-axis 
of the multi-axis accelerometer device is parallel to the Axis of Rotation of the 
turntable, and the Z-axis of the multi-axis accelerometer device is aligned 
radially, pointing along the direction of force Fl. 

To perform the calibration, the multi-axis accelerometer device is 
positioned in Orientation 1, and rotated on the turntable. In this orientation <|>(t) 
= co*t + <(>, where <j>, is the value of the angle <|> at time t=0, which is when data 
logging begins in Orientation 1. The voltages of each of the three 
accelerometers in the device are logged, while the turntable is rotating. This is 
repeated with the multi-axis accelerometer device in Orientations 2 and 3. 
Thus for Orientation 2, <|>(t) = co*t + <|> 2 , where <|> 2 is the value of the angle $ at 
time when data logging begins in Orientation 2 and similarly the equation for 
orientation 3 is <[>(t) = co*t + <J> 3 , where <J> 3 is the value of the angle 4> at time data 
logging begins for Orientation 3. 

In Orientation 1, the voltage output, denoted , of the accelerometer 
A in multi-axis accelerometer device 110 (e.g., corresponding to the 
acceleration measured by accelerometer A) is equal to the sum of the 
accelerations due to forces Fl, F2 and F3 acting on accelerometer A. The 
acceleration due to force Fl on accelerometer A is a A *g*sin(9)*A x *cos(<(>(t)) , 
where a A is the scale factor of the accelerometer A in multi-axis accelerometer 
device 110, and [\ Ay AJ is the alignment vector of accelerometer A in the 
frame of reference of the multi-axis accelerometer device 110. Similarly, the 
acceleration due to force F3 on accelerometer A is a A *g*sin(9)*Ay*sin(<|>(t)). We 
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denote the constant acceleration acting on accelerometer A (e.g., due to the 
constant force F2, and due to the constant component of gravitational 
acceleration parallel to the axis of rotation) by the constant K A V Thus, 
V A(1 =a A *g*sin(9)*{A x *cos(<t)(t)) + A/sinC^t))} + K A1 

Using similar notation for the other two accelerometers, we have 
V B(1 =a A *g*sin(9)*{B x *cos«)(t)) + B y *sin(4>(t))} + K B1 
V C(1 =a A *g*sin(B)*{C x *cos(<|)(t)) + C y *sin(<t>(t))} + 1^ A 

Similarly, for Orientation 2, shown in Figure ID, we have: 
V A , 2 =a A *g ii; sin(e)*{A x *cos((|)(t)) + A/sinCtft))} + K A>2 
VB,2=«A*g :i: sin(0) :f: {B x *cos((t)(t)) + B y *sin(4>(t))} + K B2 
V c , 2 =a A *g ;i: sin(e)*{C x *cos((l)(t)) + C y *sin(<|>(t))} + 

Also, for Orientation 3, shown in Figure IE, we have: 

V A ,3=a A *g*sin(e)*{A x ;ii cos((t)(t)) + A/sin^t))} + K A3 
V Bi 3==a A *g*sin(9)*{B x *cos(<|>(t)) + B y *sin(0(t))} + K B 3 
V c , 3 =aA Vsin(6)*{C x *cos((t)(t)) + C y *sin(<|>(t))} + 

Referring now to Figure 2, the components of the multifunction 
processor 120 in accordance with one embodiment of the present invention are 
shown. For clarity, the single input CI is shown, however it should be 
understood that the signal CI is comprised of the output signals of the multiple 
accelerometers V A V B and V c as depicted in Figure IB. 

In Orientation 1, the multifunction processor 120 first uses a low pass 
filter 201 to filter the voltage outputs (e.g., , V B1 , V C1 ) of the multiple axis 



TRMB 964 



15 



July 27, 2001 



5 of accelerometer 110. An analog to digital converter 202 converts the filtered 
signals received from low pass filter 201 into digital form. This process of low 
pass filtering and then performing the analog to digital conversion on the 
signals is referred to as sampling. The number of times a second the sampling 
is performed is referred to as the sampling frequency (f 3 ). Low pass filter 201 is 
10 designed to block all signals that have a frequency greater than a certain 

frequency, referred to as the cutoff frequency of the low pass filter, but to let 
through, or pass, all signals that have a frequency lower than the cutoff 
frequency. In the present embodiment, low pass filter 201 has a cutoff 
frequency that is lower than fJ2. 



.0 15 



The computer system 404 then logs this sampled data. This logging is 
performed for a period of time longer than that required for multiple complete 
^ revolutions of the axis of rotation. Computer system 404 then takes a 

iSi - - _ 

Discrete Fourier Transform (DFT) of each of the logged voltage outputs from 

W 20 accelerometers A, B and C. Similar low pass filtering, analog to digital 

**^| 

O conversion, data collection, and subsequent Discrete Fourier Transform is 

performed with the accelerometer in each of the three orientations depicted in 
Figures 1C, ID and IE. 



25 Consider the accelerometer A, in Orientation 1. The voltage on 

accelerometer A is given by: 

V A ,i=a A *g*sin(e)*{A x *cos((t)(t)) + A y *sin«)(t))} + K A1 
In the Fourier domain, all the energy in constant force F2 and the 
constant component of gravitational acceleration parallel to the axis of 
30 rotation, corresponding to the constant term K M is seen around zero 
frequency. Also, all the energy in the time varying forces Fl and F3, 
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corresponding to the terms a A *g*sin(9)*A x *cos(<(>(t)) and 

a A *g*sin(6)*A y *sin((|)(t)) respectively, is seen concentrated in the discrete DFT 
bins corresponding to the frequency co. This can be seen as a peak in the DFT 
at the bins corresponding to the frequency co. Since the DFT is a linear 
operation, the amplitude of the DFT at the frequency co is proportional to 
a A *g*sin(0). Where the peak value of the DFT at the frequency co is referred to 
as P M , the following equation holds true: 

P A1 = K*a A *g*sin(e)*exp(i*<[) 1 )*( \- i*Ay) 
In this equation, k is a constant of proportionality obtained during the process 
of taking the DFT of the signal V M . Similarly, the equations from the other two 
accelerometers in Orientation 1 are: 

P B1 = K*a B *g*sin(e)*exp(i*(() 1 )*( B x -i*B y ) 

P C1 = K*a c *g 5i: sin(0) ;is exp(i*(t) 1 ) :it ( C x -i*C y ) 

Using this same process the equations for the peak DFT values for 
Orientations 2 are: 

P^ 2 = K*a A *g*sin(6)*exp(i*<t> 2 )*( i*^) 
P B2 = K*a B Vsin(e)*exp(i*(^)*( B x + i*B 2 ) 
P C2 = K*cx c *g ;i: sin(e) :t: exp(i*(l) 2 )*( C x + i*C z ) 

Similarly, the equations for the peak DFT values for Orientations 3 are: 
P A3 = K*a A *g*sin(0)*exp(i :i: (|)3) :i! ( i*Ay) 
P B 3 = K*a B *g :l: sin(0)*exp(i*<l)3)*( B z + i*B y ) 
P C 3 = K*a c *g*sin(e) 5i: exp(i*<l)3) :i; ( C 2 + i*C y ) 

The computer system 404 also generates a sampled version of a sine 
wave internally, of amplitude a_nominal* g_nominal* sin(9_measured)* 
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5 cos((j)(t)). This sine wave comprises the predicted output of the accelerometer 
110. Here a_nominal is the expected scale factor of the accelerometer, 
g_nominal is the value of earth's gravitational acceleration at the point the 
accelerometer measurements are being taken, 9_measured is the measured 
tilt of the axis of rotation, and <|>(t) is a monotonically increasing function with 
10 constant derivative (co). Computer system 404 takes the DFT of this sampled 
sine wave. The energy in this sampled sine wave can also be seen in the DFT 
as a peak in the DFT bins corresponding to the frequency of the sine wave. 

J The energy in this sampled sine wave is also seen as a peak of this DFT 

15 at the bins corresponding to the frequency co. Again, since the DFT is a linear 
operation, the amplitude of the DFT at the frequency co is proportional to 
IJ a_nominal*g_nominal*sin(9_measured). Where the peak value of this DFT is 

V!! 

"I! 

referred to as P_nominal, the following equation holds true: 

Q 

Sj P_nominal = K*a_nominal*g_nominal*sin(9_measured) 

Fiji 

%| 20 where k is the same constant of proportionality found in calculating P M . 

The tilt angle 9 is measured using an accurate tilt sensor. Therefore, 
9_measured = 9. Also, g_nominal = 9.80665 m/s 2 = the approximate 
acceleration due to gravity. The errors in measuring the tilt angle and gravity 
25 are negligible in comparison to the required accuracy of the scale factor and 
alignment computation, and hence, these errors are ignored. 

Using this equation for P_nominal and the equations for the peak DFT 
values we can solve for the scale factors a A , a B> and a c and further solve for the 
30 alignment vectors [A^A^AJ, [B x B y BJ and [C^C y CJ. One method to solve 
these equations is shown in Appendix A. 
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Additionally, variations of this method, requiring fewer orientations of 
the multi-axis accelerometer device, can be implemented. For example, 
variations can be made if some information about the device is known prior to 
taking the calibration measurements, or if lower accuracy is required in the 
10 calibration. 



O 



One example of a variation is the following. If one of the sensitive axes of 
the device is known to lie in the X-Y plane in the frame of reference of the 
device, then the alignment angles and the scale factors for this sensitive can be 



- ! tii 

% 15 computed using data recorded only in orientations 1 and 2. 



Vii 
-a 

Nl 

fy|20 can be performed using only Orientation 1, using a combination of Orientation 

Nl 

fi 1 and Orientation 2, or a combination of all three Orientations. The errors in 

1-4 



Another example of a variation, is the case in which one of the sensitive 
axes lies nominally in the X direction of the frame of reference of the device, and 
the accuracy requirements on the sensor are lower. Then, a similar procedure 



the scale factor and alignment angles of the sensitive axis become smaller, in 
this example, as more orientations are used for the calibration. 

25 Both of these variations can be solved as approximations to, or special 

cases of, the general solutions described in Appendix A. 



Referring now to Figure 3, a flowchart showing the steps of a process 
300 in accordance with one embodiment of the present invention is shown. As 
30 depicted in Figure 3, process 300 shows the operating steps of the calibration 
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system (e.g., system 100 of Figure 1A) calibrating a multiple axis multi-axis 
accelerometer device (e.g., multi-axis accelerometer device 110). 

Process 300 begins in step 301, in which a tilt angle of the turntable 
mechanism, on which the calibration of the accelerometer device is to be 
performed, with respect to the local gravity vector is measured. As described 
above, this tilt angle (e.g., 9 ) describes a difference between an axis of rotation 
of the turntable and the local gravity vector (e.g., the vertical axis). In step 

302, the multi-axis accelerometer device to be calibrated is mounted on the 
turntable mechanism. The mounting is performed such that the Z-axis of the 
multi-axis accelerometer device is pointing along the axis of rotation. In step 

303, the turntable is spun around the axis of rotation at an angular velocity co. 
As described above, the rotation gives rise to three forces acting on the 
accelerometer; a centripetal force F2 of constant magnitude acting on the 
accelerometer, and time varying forces Fl and F3 due to components of the 
local gravity vector. The constant component of gravitational acceleration, 
acting along the axis of rotation, also acts on the accelerometer, and is not 
shown in the figure. A component of Fl and F3 is simultaneously experienced 
by each accelerometer (e.g., A, B and C) of multi-axis accelerometer device 
110. In step 304, the V B and V c outputs of the accelerometer device are 
filtered using a low pass filter. In step 305, the filtered outputs of the 
accelerometer device are processed through an analog to digital converter 
yielding respective sampled outputs of the accelerometer. In step 306, the 
sampled outputs of the accelerometer device are logged as the accelerometers 
A, B and C of the multi-axis accelerometer device experience the time varying 
components Fl and F3. 
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Referring still to process 300 of Figure 3, in step 307, the multi-axis 
accelerometer device is detached from the turntable and remounted such that 
the Y axis of the device is pointing along the axis of rotation. In step 308, steps 
303 to 306 are repeated with the multi-axis accelerometer device in this new 
orientation. 

In step 309, the multi-axis accelerometer device is detached from the 
turntable and remounted such that the X-axis of the multi-axis accelerometer 
device is pointing along the axis of rotation. In step 310, steps 303 to 306 are 
repeated with the multi-axis accelerometer device in this new orientation. 

In step 311, a predicted output of an ideal accelerometer parallel to the 
plane of rotation of the turntable is generated. In step 312, discrete Fourier 
transforms are performed on all of the logged accelerometer outputs and the 
predicted output of the ideal accelerometer. Hence, in step 313, the discrete 
Fourier transforms of each the logged outputs along with the predicted output 
of the ideal accelerometer are used together to determine the scale factors and 
alignment vectors of each of the accelerometers of the multi-axis 
accelerometer device. 

In so doing, the turntable mechanism of the present invention 
accurately measures and determines the scale factors and alignment vectors 
of the multiple axis accelerometer device without relying on any time varying 
control of a standard device (e.g., stepper motor, etc.) to impart variable 
acceleration to the accelerometer. Furthermore, it determines these scale 
factors and alignment vectors without the knowledge of the distance to any 
measurement point of the multi-axis accelerometer device from any other 
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point. 



Thus, the present invention provides a solution that accurately 
measures and determines the scale factor and alignment vector of each 
accelerometer in a multiple axis multi-axis accelerometer device 
simultaneously. The present invention provides a solution that can calibrate a 
device having multiple sensitive accelerometer axes in a single calibration 
process. The present invention provides a solution that calibrates multi-axis 
accelerometer devices without introducing unnecessary sources of error. The 
solution of the present invention is precise and avoids reliance on standard 
devices, which can introduce error into the calibration process. The solution of 
the present invention does not rely on any time varying control of a standard 
device to impart variable acceleration. Furthermore the solution of the 
present invention does not rely on a precise measurement of the distance to 
any measurement point multi-axis accelerometer device from any other point. 

Computer System Environment 

Figure 4 shows a diagram of the basic components of computer system 404. 
Within the above discussions of the present invention, certain processes and steps 
are discussed that are realized, in one embodiment, as a series of instructions (e.g., 
software program) that reside within computer readable memory units of 
computer system 404 and executed by the processor(s) of system 404. When 
executed, the instructions cause the computer system 404 to perform specific 
actions and exhibit specific behavior, which was described in detail above. 

In general, computer system 404, used by the present invention, comprises 
an address/data bus 451 for communicating information, one or more central 
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processors 453 coupled with bus 451 for processing information and instructions, a 
computer readable volatile memory unit 452 (e.g., random access memory, static 
RAM, dynamic RAM, etc.) coupled with bus 451 for storing information and 
instructions for central processor(s) 453, a computer readable non-volatile 
memory unit 454 (e.g., read only memory, programmable ROM, flash memory, 
EPROM, EEPROM, etc.) coupled with bus 451 for storing static information and 
instructions for central processor(s) 453. Computer system 404 interfaces with 
the other components of system 100 via system interface 460. System 404 can 
optionally include a mass storage computer readable data storage device 455 , 
such as a magnetic or optical disk and disk drive coupled with bus 451, for storing 
information and instructions, a display 456 for displaying information to the 
computer user, and an input output device 457 including, for example, 
alphanumeric and function keys for communicating information and command 
selections, cursor control inputs, command selections, etc. 

The foregoing descriptions of specific embodiments of the present 
invention have been presented for purposes of illustration and description. 
They are not intended to be exhaustive or to limit the invention to the precise 
forms disclosed, and obviously many modifications and variations are possible 
in light of the above teaching. The embodiments were chosen and described in 
order best to explain the principles of the invention and its practical 
application, thereby to enable others skilled in the art best to utilize the 
invention and various embodiments with various modifications as are suited to 
the particular use contemplated. It is intended that the scope of the invention 
be defined by the Claims appended hereto and their equivalents. 



TRMB 964 



23 



July 27, 2001 



